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Abstract 

We derive the first order canonical formulation of cosmological per- 
turbation theory in a Universe filled by a few scalar fields. This theory 
is quantized via well-defined Hamiltonian path integral. The propaga- 
tor which describes the evolution of the initial (for instance, vacuum) 
state, is calculated. 

1 Introduction 

The quantization of linearized matter and metric perturbations in the Fried- 
mann Universe has become an especially interesting topic in connection with 
inflationary scenarios of the Universe evolution. 

Actually, if the Universe went through an inflationary stage, then the 
initial quantum fluctuations can explain the observable large-scale structure 
of it (see, for instance, Q). 

One of the simplest models for inflation is the chaotic inflation when 
the Universe is filled by homogeneous scalar fields [Q]. In this model, small 
inhomogeneities in the distribution of scalar fields create the metric pertur- 
bations. In addition there can be fluctuations of the metric which are not 
due to inhomogeneities of the matter (gravitational waves). 

The purpose of this paper is to present a selfconsistent quantum theory 
of these perturbations starting from first principles. 

It was shown in [fjj that the Lagrangean for the scalar cosmological per- 
turbations derived via expansion of the Einstein-Hilbert action can be ex- 
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pressed entirely in terms of one gauge-invariant variable which describes the 
collective degree of freedom of the metric and matter perturbations. This 
reduced Lagrangean was then the starting point for the quantization. 

There are several disadvantages in the above mentioned semi-quantum 
approach. First, it does not permit us to reveal the connection between 
fundamental quantum theory of gravity and quantization of small metric 
perturbations. In particular, the role of gauge group (diffeomorphisms) is 
not clear there. Second, to reduce the action to the function of only one 
variable, we need to use some of the Einstein equations. 

It is much more easy to clarify the role of these equations in the canonical 
quantum theory, where they are just algebraic constraints. Further, the 
calculations in Lagrangean theory are not very straightforward and it is a 
tricky point to generalize them to include in the consideration several scalar 
fields i|. 

And finally, the Lagrangean theory is not very suited to relate the func- 
tional approach to quantization of perturbations with "fundamental" canon- 
ical or path integral quantum gravity where the question about "boundary 
conditions" for a quantized Universe is well posed || 0. 

Our method in this paper is to expand the canonical ADM action for 
gravity and N scalar fields up to second order in the perturbation variables. 
Thus the first order formulation of the classical theory of cosmological per- 
turbations is derived. Then we quantize this theory via path integral in the 
Hamiltonian formalism where the measure is well defined. As a result the 
integration over the (infinite) volume of group (diffeomorphisms) is factor- 
ized explicitly and the problem is reduced to the quantization of a set of 
gauge invariant fields with time-dependent masses. The propagators and 
evolution of initial vacuum state are calculated. The details of calculations 
will be presented in a forecoming paper Q. 

We use the units in which c = Ti = 16nG = 1 and adopt MTW conven- 
tions [|. 



2 Action 



If we write the metric in the ADM form [10|, 



ds z 



■{M 2 - MiM^dt 2 + 2M i dx i dt + jijdx^x^ 
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then the first order action for gravity and iV-scalar fields if a (A = 1, • • • , N) 
with potentials V A {y) is g, [□]: 

S = [ if'.nlti - 1 ^ ,, + £ n A <p A - N a H a ), (2) 
J A 

Here, the dot means the derivative with respect to coordinate time t, ti 13 
and ir A are the momenta conjugated to %j and ifA- The lapse Af(x,t) and 
the shift M l (x,t) functions play the role of Lagrangean multipliers. 
Correspondingly, the four constraints are 

Ho = G ijkl ir i iir kl -^WR+ (3) 



.4 



2^7 2 



A 



where 



G ijk i = — — (7ife7j7 + liiljk - Hjlki) , 7 = det(7y) (5) 

7 

and the bar | denotes covariant derivatives with respect to the three-metric 

lij- 

3 Background 

In an isotropic, flatQ Universe, the coordinate system can be choosen in such 
a manner that 

¥A = <PA(t), ti> = n A (t), 

lij = a 2 (t)Sij, tt 1 J> = <%7r a /6a (6) 
M = M(t), M* = 0. 

Then the action @ reduces to 

5 = | dtd 3 x (vr a a + rf<PA ~ AAW ) , (7) 



x The results can easily be generalized for closed and open Universes. 
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where now the superhamiltonian 

\2 /UA\2 



is constrained to vanish. Using Hamilton eqs. we can express the momenta 

V 



7r a and tt a in terms of the "velocities" a and ip A : 



aa a a 3 

*a = -12 jf, 7T A = jjtfA. (9) 

From the vanishing Hamiltonian (^) and other Hamilton eqs. with taking 
into account @, one obtains the the following eqs. for the background 
Friedmann Universe: 

Cp A + (3H - !t)<p A + M 2 V A = 0, (10) 



6 ^ 2 = EQ(^) 2 +aaV^)), (ii) 



where H = a/a coincides with the Hubble parameter only if we choose the 
gauge M{t) = 1. Later, these eqs. will heavily be used to simplify the action 
for cosmological perturbations. 



4 Perturbations 



To study the behaviour of perturbations in a flat Universe we consider small 
deviations from the homogeneous solutions @, such that for the 3-scalars 
one has 



<p A (x,t) = <p A (t)+5tp A (x,t) 

7T A (x,t) = 7T A (t) + 57T A (x,t) (12) 

M(x,t) = M(t) + 5M{x,t) 

It is convenient to write the perturbations in 3-vector 5M l and 3-tensor 6jij 
and as a sum of scalar (s), vector (v) and tensor (t) peaces which don't 
interfere in the linear approximation: 

SM i = (S) B' 1 + {v) D i 

5^ = a 2 [-2 (s V<% + 2 w E 4j + MVij + (v) Vj,i + 

5n ij = J-^D^' + (s2) vr' ij + <V J ' + 'V'* + V] (13) 
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where the comma means usual derivative with respect to corresponding spa- 
tial coordinates. The vectors D l ,Vi,it l and tensors (t) hij, (t) ir l:) satisfy the 
following conditions: 2 

(^. = 0, (V, = 0=^ 1 1 

We will raise and lower indices by the unit tensor 5{j. 

Expanding the action (Q) in perturbations, we find that the first order 
terms vanish on the eqs. of motion for the background model whereas the 
second order terms give the action for the perturbations. 

This action consists of three decoupled terms which correspond to scalar, 
vector and tensor (gravitational waves) perturbations: 

5 2 S = «S(<VV<V^, (sl) ^ (s2 M 

The explicit form of the action for scalar perturbations up to total deriva- 
tives is 



S = J dtd 3 x U^ + npF + Y, 5 rf 5l PA - (a) A 0(s) C - (s) A x {s) d - (s) H 

(16) 



.4 

where the two constraints linear in momenta are 



"Co = -^-^{H 2 + H-H^){^-F)-Aa^ 

and the Hamiltonian for scalar perturbations (s) W takes the form 

{B) n = ^3(34 + 27r F 7r^)-F(7r^V + ^ J F + 4vr F F) 

4a 3 Af 
+ -jjr (H 2 — H + H—)(3^ 2 - 2^F + F 2 ) + 2NaAipip 



2 Here and in what follows, we omit repeating the superscripts (s) and (v) whenever it 
is clear which of the decoupled peaces is considerd. 
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+ £ Up A &4 - Ma 3 V$5ip A m — F) + ^ («5<) 2 

-^AVW^^W'), (18) 
and we introduced the new independent variables^ 

F = AF, 

up = 2(( sl >7r + A< s2 >7r), = -2(3< sl 'vr + A^vr), (19) 

«A° = 4A/", ^A 1 = AS. 

For the vector perturbations one gets: 

c«>5 = y rftd 3 !^^ - W A* W Q - (20) 

where 

W W = -j r (H 2 -H + H J -^)V i AV l -2H7r i v V i , (21) 

and we introduced instead of (l,) 7r* and D l the independent variables 

Trir = -2A (l V, w A* =£) i + ^ ( V (22) 

a 6 

Note that in this case we also have only two independent constraints, since 

M X\i = , = 0. (23) 

For the tensor perturbations there are no constraints and the action takes 
the form 

W S = f dtd 3 x(ir%hij - (t) n), (24) 



where irV = ^tt and 



{t) n - —TT ij TT ij -AHTr ij h-- 



a" ,„ -Jv\, , Ma 



4- jj(H 2 -H + H-jfihijha - '-fAhijhij. (25) 



3 Here and throughout the paper, A denotes the flat three-dimensional Laplacian, e.i.. 
AE = E.u. 
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Let us mention an interesting point. Starting with a superhamiltonian and 
fixing the background, we end up with a pure Hamiltonian for perturba- 
tions which is responsible for their dynamics with respect to background 
time. There are still four independent constraints (two for scalar perturba- 
tions plus two for vector ones), as it should be. These constraints are linear 
in momenta and generate themselves the gauge transformations which cor- 
respond to diffeomorphisms fl2|j . 



5 Gauge-invariant variables 

The constraints (s) Co, <s) Ci, M Cj and the Hamiltonian 

Hit) = J d 3 x( (s) H+ <V) H+ {t) H) 

form a closed algebra with respect to Poisson brackets {..,..} defined in the 
standart manner for our canonically conjugated variables 5(p a , Stt^; ip, tt^; F, 
7rp; Vi, n % v ;hij,ir%; that is 

lc a ,Cp\ =t^C 7 and ^H(t),C a \=4C + ^. (26) 

For instance, for the scalar constraints (s) Co, (s) Ci we have 

C (x,t),C (y,t)} = = {c 1 (x,t),C 1 (y,t)y 

Cb(M),} = ^A x C 1 (x,t) + dCo ^ t) , (27) 



H(t),d(x,t) 



dCi(x,t) 
di ' 



The action (16) is invariant with respect to the transformations generated 
by the constraints (|l7]) , if we simultaneously transform the Lagrangean mul- 
tipliers (a) A": 

% = U,t a c a } 



where £ = £ a are the parameters of transformation and q can be any of the 

V 



canonical variables 5ip a , 5ir A ij) , etc. 



4 The appearance of time derivatives on the right-hand-side is due to the explicit time 
dependence of the constraints. 
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Since the constraints are linear in momenta, the transformations (|2 



correspond to diffeomorphisms [Oj. From (28) and ( |16D with taking into 



account pq), (|27|), it is easy to get the following transformation laws: 

^A 1 = i l - MAS, /a 2 , <%A° = i°, (29) 
<%F = I 1 , etc 

Using ([29| ) one can easily construct gauge invariant variables Q for which 
5^Q = 0. For instance, the most important of them are: 

v A = a{5ip A + (30) 
A<D = I A A°-i|-[^(A 1 -F)]. (31) 

As we will see, the first variable naturally appears when we quantize the 
theory, whereas the second one has clear physical interpretion, namely, <3? 
corresponds to the relativistic Newtonian potential jlj. 

The gauge invariant variables for the vector perturbations can be con- 
structed in the same manner. The tensors h+j and 7Tif are gauge invariant 
by themselves since there are no constraints which contain these variables. 
For more details, refer to jq]. 



6 Path integral 

To quantize the perturbations we write the path integral (h = 1) 

K = J* V^S+^S+US) (32) 

which is just the propagator from some initial (i) to final (/) state. The 
diffeomorphism invariant measure T>\x can be taken as 

Vfi = V(5j tj )V(5ir ij )V(5N a ) = J P w fW M fJ,T> (t) // (33) 

where J is an irrelevant constant Jacobian which can be absorbed in nor- 
malization and 

V (v) fi = VV l Vir\rV (v) \\ (34) 
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are the measure, correspondingly, for scalar, vector and tensor perturba- 
tions. For T>(8(p A ),T>(STT A ),T>ip,T>Tr^ etc, we take the usual Liouville mea- 
sure. Then the propagator K can be represented as a product of propagators 
for different types of perturbations. Let us calculate them separately. 
For the scalar perturbations: 



K 



V 



fie 



(35) 



where is given by (16) and P (s) fi is defined in (34). 

Let us integrate at first over the Lagrangean multipliers ("'A /"' A . Then 
we get the delta functions of constraints (s) C a , and since they are linear in 
momenta they permit us to integrate easaly over tt^,7Tf. Then, chang- 
ing the variables (in particular, introducing the gauge invariant variables 
va, see (|30|), instead of S(p A ), with help of the eqs. of motion for back- 
ground (||)7 ( |io|) and (|Tl| ) and after integration over the momenta we obtain 
the following result^] 



{S) K = M 



Y[Vv A e 

A 



S(v A ) 



t>v A =h 



dv{ 



where 



M oc / VipVF 



(36) 



(37) 



corresponds to the explicitly factorized infinite volume of the gauge group 
and can be absorbed in normalization. The action S(va) depends only on 
the gauge invariant variable va, (see (|30|) ) : 



S(v A ) = \f d 3 xdJj2(v'i - vjvj) + Vab 

Z J \ A /ID 



v A v B 



(38) 



A A,B 
Here the prime means derivative with respect to the new time parameter 

Mdt 



V 



and the function of time £Iab is 



n 



AB 



^ + 2(_ 
za H 



H) 



+ 



ZA 



a^A 
H ' 



ZA 



El 

H 



(39) 



2H< 



Sab + -^(Hzazb)' 



(40) 



3 We also skip from the action a lot of total derivatives. 
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We see that the metric perturbations ip, F do not have their own dynam- 
ical degrees of freedom. They are entirely due to the perturbations in the 
matter distribution. In this sense the quantization of scalar cosmological 
perturbations is just the same as quantization of matter itself with taking 
into account the gravitational field created by this matter. 

The vector perturbations do not possess any dynamics. This can eas- 
ily be seen if we integrate in ^S, at first over the Lagrangean multiplier ^X 1 
and then over the momenta ir l v . As a result we find that the action {v) S{Vi) 
vanishes (up to the total derivatives). The integral over Vi corresponds to 
the integration over the gauge volume, like in ([37]). 

If the Universe would be filled by some matter which has vector degrees 
of freedom (like, for instance, by a perfect fluid) then the quantization of 
vector perturbations would become nontrivial. 

For tensor perturbations the calculations are very simple. Integrating 
over the momenta tt 1 ^ one arrives at: 

®K = J VliijP-^ 1 (t)s = J Ve ije l (t)s ^ , (41) 

where 

1 f 3 / / a " 

S — — J dfjd x[e^je^j &ij y m&ij,m H i (^^) 

and ejj = ahij/V2 is the the gauge invariant variable. 

Thus the quantization of the cosmological perturbations in a Friedmann 
Universe has been reduced to the quantization of a set of gauge invariant 
fields with time-dependent masses in a flat space-time. 



7 Propagator 

As an example, we calculate the propagator and the evolution of an initial 
vacuum state for scalar perturbations in the Schrodinger representation. To 
simplify the consideration we assume that there is only one scalar field in 
the Universe, that is va = v and VLab = = z" / z. 
Expanding v(x,rj) in Fourier modes as 



n) = ^(j (2^)3/2 ^) Sm (^) + / (2^)3/2 ^) COS (^)) ( 43 ) 



v(x. 



( 2vr )3/2 "VU^K J (27r)3 /2 

fe 3 >0 fc 3 <0 



where the Fourier coefficients vt are real, one gets 

S(v) = \J d 3 kd V (v' k 2 - k 2 v 2 k + 0(7?)^. (44) 
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Thus, we can write the propagator from some initial (i) to final (/) field 
configuration as 

K(f\i) oc Y[K k {v f ,rif\vi,rH) t (45) 

k 

where| 

K k (v f ,rj f \vi,r]i) = J Vv k expjj J dr](v' k 2 - k 2 v 2 k + ftfaK)} ( 46 ) 

is the propagator for the harmonic oscillator with the time-dependent fre- 
quency 

Lo 2 k ( V ) = k 2 -tt(r 1 ). (47) 
The result for the path integral (|4FJ) is well known [13, [T|] : 



^*«> = fip;^ M ' (48) 

where S^/ is the action for the classical trajectory satisfying the equation 

v'l + u? k v k = (49) 

and the boundary conditions v k (rji) = Vi,v k (r]f) = Vf. 

The classical action S c i can be expressed as a combination of two funda- 
mental solutions C{rf),D{rj) of the eq. fl49]); 



v k (v) = ViC(ri)+VfD fa) (50) 

Here the functions C and D depend on k, but not on Vi, Vt and they satisfy 
the boundary conditions 

C(r 7i ) = l, C(r?/)=0, 
£>fa) = 0, J D(r ?/ ) = l. 



In terms of C and D, the explicit expression for the propagator ( |46[) is given 
by: 

1 X 1 / 2 



(52) 



6 To simplify the notation we also skip the index in some of the formulae whenever it 
is clear which of functions depend on k. 
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where 



f(Vf,Vi) 



d 

and Ca ee —C(r]i) ,etc 



(53) 



do not depend of Vi,Vf. 

Given the propagator, we can calculate the evolution of any initial state 
functional ^>(vi,rji). For instance, if this initial state can be represented in 
the form 

*( Vi(x, r]i), rj ) ocJJ* fe (wj, rji) , (54) 



then as a result of evolution we find that 



<S[v f (x,r)f),r]) oc Y[y k (v f ,r)f) 

k 



where 



+00 



Vk(vf,T)f)= / dviK k (vf,rjf\vi,rn)^> k (vi,rji). 



(55) 



(56) 



Thus, there will be no mixing of the Fourier components with different k in 
the linear approximation as it should be. 

If we choose, for instance, a normalized gaussian state with dispersion 
ctq (which can depend on k) as initial wavefunction, 



vM^H^r^expl-^} 



2a 2 J ' 



then the result for the integral (56) is 

^( V/ ,r ?/ )=e^(vrS 2 )- 1 /4 exp |_^2 ( ^__ zA) 

where S, A and 9 depend on i]i,rjf, k and ctq: 



S 2 EE 



0"n 



1 + oiC? 1 



W - c f ) 2 



(57) 



(58) 



(59) 



COt 6 EE (TqC[. 



Thus, the wave function at time r)f describes still a gaussian state, but now 
with a modified complex dispersion. 

Since all information about the system is contained in the state func- 
tion (|5"8|), it is standard to compute expection values (i>kVk), etc. 
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If we want to take as initial state the vacuum one, we need to specify the 
dispersion oq entering (|57|) . As it is well known there is an ambiguity in the 
definition of vacuum in the presence of external fields . In our case this 



ambiguity is due the time dependence of the frequency tuk in (47)- However 
for any definition of vacuum we should have 

al — ► -r for k 2 > U(k) (60) 

Ki 

Later we shall see that this condition is sufficient to get unambiguous pre- 
dictions for the spectrum of fluctuations in the most interesting range of 
scales in inflationary model. 

The final step is to calculate the two-point correlation function of the 
gauge invariant relativistic potential $>(x,rj), (see (|3l"|)). 



8 Correlation function and power spectrum 

For the purpose of comparing with observation we need the two-point cor- 
relation function £(r)[] of the gauge invariant metric fluctuation &(x,r]): 



£(r) = rj)\$(x)$(x + r)\t>(v, 17)) (61) 

-N 2 (62) 



dk s\n{kr) 2 



kr 







where ^(v,rj) is the state functional and the power spectrum \5k\ is the 
measure of the amplitude of metric fluctuations on a comoving scale ~ 1/k. 
The relation between <E>(x, 77) and v(x,r/), and, correspondingly, between the 
quantum operators 4? and v, follows from Hamilton eqs., (see 0): 

A$ = ^-(n v --v) (63) 

where tt v is the momentum conjugated to v. 

In the Schrodinger representation, working with separate Fourier modes, 
one has for the appropriate momentum operator % = —id/dv^. 



Now using the definitions (|6l|), taking into account (62) and (|63| ) and 
calculating the corresponding correlation functions 

(VkVk) , (TTfcTTfc) , (VkKk + TTfcTJfc) , (64) 



7 In the case under consideration the correlation function just depends on r — \r\. 
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for the state fl58D, we get the following result for the power spectrum \5k\- 



Thus, given an initial gaussian state at i] = r/i, specified by the dispersion 
(Jo(k), we derived the power spectrum \5k\ at an arbitrary later moment of 
time r] = r]f. This spectrum can be expressed entirely in terms of the 
two functions C(rj) and D{rj) which satisfy the eq. ( |49|) with the boundary 
conditions (|5l|). 

The functions X and A, in terms of C and D, are given in (|59|). For an 
initial "vacuum" state the asymptotic dependence of <Jo(k) on k at k 2 3> Q 
is specified in (|60|). 

As an example, let us consider chaotic inflation in the model with a 
massive scalar field: V((p) = ^m 2 ip 2 . 

In this case from (|65| ) one obtains the following spectrum for the metric 
fluctuation after the end of inflation in the most interesting range of galactic 
scales: 

|4| ~ (V3/107r)mln(A p/l /A 7 ), (66) 

where X p h = a(rj)/k is the physical wavelength of the perturbation and A 7 
is some characteristic wavelength of the cosmic background radiation. This 
result coincides with the one obtained by other methods (see Q). The in- 
flation amplifies the initial vacuum fluctuations enough to explain the large- 
scale structure of the Universe only if m ~ 10 15 GeV. The other interesting 
applications are considered in a forecoming paper ]|]. 

9 Discussion 

Starting with the Hamiltonian (ADM) formulation of General Relativity, 
we deduced the Hamiltonian theory of cosmological perturbations in the 
Friedmann Universe. 

To derive the action for linearized metric and matter perturbations in a 
first order Hamiltonian formalism, we expanded the ADM action for gravity 
and the action for the matter (scalar fields) up to second order in pertur- 
bations. The concrete calculations have been done for the case when the 
background Friedmann Universe has zero spatial curvature. This is a good 



47T 

(2tt) e 



1 



8ira k 



(65) 
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approximation if, afterwards, we want to quantize only the perturbations 
assuming that the background is classical. 

The generalization of the developed formalism to the closed Universe is 
straightforward and will be presented in the future publication || . This case 
is especially interesting if we also want to quantize the Universe as a whole. 
However if the zero mode perturbation which correspond to the Friedmann 
background is in a quasiclassical region, then only the quantization of inho- 
mogeneities becomes interesting. 

We have further studied the diffeomorphism transformations in Hamilto- 
nian formalism and constructed explicitly gauge (diffeomorphism) invariant 
variables for perturbations. 

Then, the first order Hamiltonian theory was quantized via well defined 
Hamiltonian path integral. The volume of the gauge group (diffeomorphism) 
was factorized explicitly and the problem was reduced to the quantization 
of a set of gauge invariant fields and, finally, to the quantization of harmonic 
oscillators with time dependent frequences. 

We calculated the propagator which describes the evolution of the quan- 
tum state and then found a closed expression for the power spectrum of 
fluctuations for initial gaussian, for instance, vacuum state. 

All the calculations have been done in the Schrodinger picture which 
is very convenient to study, for example, the decoherence problem for the 
cosmological perturbations ju], [D]]. The results were finally applied to the 
concrete model of chaotic inflation. 
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